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Summary 

A water tunnel investigation was conducted to 
demonstrate the capabilities of a laser-based in- 
strument that can measure velocity and fluores- 
cence intensity simultaneously. Fluorescence inten- 
sity of an excited fluorescent dye is directly related 
to concentration level and was used to indicate the 
extent of mixing in flow. This instrument is a 
three-dimensional laser Doppler velocimeter (LDV) 
in combination with a fluorometer for measuring flu- 
orescence intensity variations. This capability allows 
simultaneous flow measurements of the three orthog- 
onal velocity components and mixing within the same 
region. Two different flows which were generated by 
two models were studied: a generic nonaxisymmetric 
nozzle propulsion simulation model with an auxiliary 
internal water source that generated a jet flow and 
an axisymmetric forebody model with a circular sec- 
tor strake that generated a vortex flow. The off-body 
flow fields around these models were investigated in 
the Langley 16- by 24-Inch Water Tunnel. The ex- 
perimental results were used to calculate 17 quan- 
tities that included mean and fluctuating velocities, 
Reynolds stresses, mean and fluctuating dye fluores- 
cence intensities (proportional to concentration), and 
fluctuating velocity and dye concentration correla- 
tions. An uncertainty analysis was performed to es- 
tablish confidence levels in the experimental results. 
In general, uncertainties in mean velocities varied be- 
tween 1 and 7 percent of free-stream velocity; un- 
certainties in fluctuating velocities varied between 1 
and 5 percent of reference values. The results show 
characteristics that are unique to each type of flow. 

Introduction 

To measure as many characteristics as possible 
of the off-body flow around aerodynamic shapes and 
within internal flow fields, instruments were devel- 
oped that provide flow information previously un- 
available. Mechanical probes often adversely affect 
or disturb the flows they are measuring. Therefore, 
emphasis was placed on the use of nonintrusive opti- 
cal measurement techniques which can provide useful 
qualitative and quantitative details of many types of 
flow. In particular, the development of the laser gave 
the researcher a tool for making accurate determi- 
nations of many flow characteristics without signif- 
icantly disturbing the flow field. Volumes of litera- 
ture have been written which describe techniques for 
nonintrusive optical flow measurements. Reference 1 
provides a thorough overview of general laser-based 
flow measurement concepts. 

Laser Doppler velocimetry has been used to inves- 
tigate flow fields since the mid-1960’s and is described 


extensively in the literature. (See refs. 2 and 3.) 
Fluorescence intensity measurement is a well-known 
spectroscopic technique and fluorometry concepts are 
described in reference 4. 

The nonintrusive technique used in the experi- 
ments discussed in this paper was a combined laser 
Doppler velocimetry (LDV) measurement of veloc- 
ity and laser-induced fluorescence intensity measure- 
ment of dye concentration. A system was developed 
for the Langley 16- by 24-Inch Water Tunnel and was 
first described in reference 5. The three orthogonal 
velocity components of flow were measured by the 
LDV section of the instrument; the scalar fluores- 
cence intensity was measured simultaneously by the 
fluorometry section at the same location as the LDV 
measurements and characterized mixing in the flow. 

The simultaneous measurement of three velocity 
components and the level of dye concentration in the 
sample volume is a unique capability that allows ex- 
amination of flows that have both a high degree of 
three-dimensionality and significant mixing between 
different regions of the flow field. By examination of 
the entire three-dimensional characterization of the 
flow, the predominant time-averaged mechanisms in 
mixing as well as the three-dimensional velocity flow 
field can be determined. This had not previously 
been accomplished. To demonstrate the capability 
of the instrument, two very different types of flows 
were examined: a turbulent jet flow and a vortex 
flow. Reference 5 reported on the design and oper- 
ation of the hardware and software systems as well 
as provided some limited test results. This paper 
presents detailed experimental results from measure- 
ments of each type of flow and provides an analysis 
of estimated experimental uncertainties. 

Symbols 

A area integral of vorticity 

a transformation matrix element 

B bias error 

b absorption path length, in. 

(see eq. (2)) 

bo jet half-height, in. 

C integration circuit for 

circulation 

c transformation matrix element; 

concentration of fluorescent 
dye, mg/ml 

D e - 2 laser beam diameter at e~ 2 of 

peak intensity levels, mm 


d 

distance of lens from tunnel 
wall, in. 

df 

sample volume fringe spacing, 
A 

0 . > 

2 sin ft 

dm 

sample volume diameter, 
4A / 

7rD e -2 cos ft ’ ^ 

dS 

surface element area, in 2 

e 

unit vector in LDV coordinate 
system 

F 

fluorescence; fluorescence 
intensity, V; focal length, in. 

f 

transmitting lens focal 
length, mm 

1 0 

incident radiant power, W 

i, j, k 

unit vectors in tunnel coordi- 
nate system 

K 

constant, (see eq. (3)) 

lm 

sample volume length, 
4A / 

7TjD e -2 shift’ 

N 

number of samples 

n 

index of refraction of media 

P 

precision limit, P = £ 95 S 

r 

result in uncertainty analysis; 
radius, in. 

S 

precision index 

*95 

statistical t distribution value 
giving 95 -percent confidence 
level 

u 

uncertainty; velocity, in/sec 

u 

axial velocity, in/sec 

u c 

jet exit velocity, in/sec 

u m 

jet centerline axial 
velocity, in/sec 

u n 

velocity perpendicular to 
fringes, in/sec 

u , v, w 

mean velocity components, 
in/sec 

u', v 1 , w' 

velocity fluctuations, in/sec 


ti /2 , x/t/ 2 , x/u / 2 rms velocity fluctuations, 



in/sec 

u f c! , u'c', 

covariance of velocity and 
concentration fluctuations, 
in.-V/sec 

v!v!, v f v f , w f w f 

Reynolds normal stresses, 
in 2 /sec 2 

uV, v f w f , w f u f 

Reynolds shear stresses, 
in 2 /sec 2 

V 

lateral velocity, in/sec 

W 

lateral depth in test section, 
in.; nozzle width, in. 

U) 

vertical velocity, in/sec 

X 

variable in uncertainty analysis 

X , y , z 

Cartesian coordinates, in. 

a 

angle of attack, deg 

r 

circulation, in 2 /sec 

A 

implies uncertainty when pre- 
ceding quantities in uncer- 
tainty analysis 


Kronecker delta 

£ 

absorptivity of fluorescent dye, 
in -1 (see eq. (2)) 

6 

angle in LDV coordinate 
system, deg 

ft 

laser beam half-angle, deg 

A 

wavelength of light, nm 

V 

frequency of laser light, Hz 

A 1/ 

shift or Doppler frequency, Hz 

P 

correlation coefficient; fluid 
density, lbf-sec 2 /ft 4 

0 

standard deviation 

$ 

quantum efficiency of fluores- 
cent dye 

0 

strake azimuth angle (0° at 
top center), deg 

u; 

vorticity, sec -1 

Subscripts: 


Bij 

correlation coefficient for bias 

errors 

cell 

quantity for measurement grid 


cell 


2 



i, 3 

elemental component 

ij 

elemental indices 

3 

jet 

LDV 

laser Doppler velocimetry 
coordinate system 

max 

maximum 

Pij 

correlation coefficient for 
precision errors 

pk 

peak or maximum 

r 

result 

rms 

root mean square 

RSS 

root sum square 

t 

tangential 

v, 9, b 

components in LDV coordinate 
system 

X 

streamwise component 

oo 

free-stream condition 

Superscripts: 


/ 

fluctuating quantity 

- 

mean quantity 

A 

unit vector 

Abbreviations: 


AD 

axisymmetric decay region of 
jet 

A/D 

analog-to-digital 

ANSI 

American National Standards 
Institute 

ASME 

American Society of Mechani- 
cal Engineers 

CD 

characteristic decay region of 
jet 

LDV 

laser Doppler velocimetry 

LSB 

least significant bit 

PC 

potential core region of jet 

PMT 

photomultiplier tube 

rms 

root mean square 

RSS 

root sum square 

SV 

sample volume 


Experimental System Description and 
Test Technique 

Measured and Calculated Quantities 

Four quantities are measured in the instrument 
coordinate system: A v UiV , A v u ,g, A v u ^ and F\ 

which are the three frequency components that yield 
three velocity components and the fluorescence in- 
tensity in the LDV coordinate system, respectively. 
From these untransformed measurements, the follow- 
ing 17 quantities are calculated: u, u, w; Vv^, 

; u f u\ v*v\ w f w f ; uV, v f w f , w f u f ; v!d , v f d , w f d\ 
and c, d . 

Instrumentation System 

Velocimeter . The use of two laser beams for 
velocity measurement is illustrated by figure 1. The 
planar wave fronts of the two intersecting beams form 
a set of interference fringes that are parallel to the 
line bisecting the angle between the two intersecting 
beams. Velocities are measured in the sample volume 
(SV) formed at the intersection of the two beams. 
The Doppler effect occurs because of the relative 
motion between the laser-scattering particle and both 
the laser source and receiver. Equation (1) can be 
derived from the principles of Doppler theory (refs. 2 
and 3) and relates the Doppler shift frequency Ai/ to 
the velocity of the particle by 


The Doppler signal is characterized by a varying 
amplitude and a constant frequency. The locus of 
the signal half-amplitudes forms a bump called the 
pedestal. This pedestal is the result of the Gaussian 
power distribution across the laser beam that peaks 
at the beam center and is generally filtered out early 
in the signal-processing procedure. 

The signal will produce a frequency that can be 
directly related to velocity. However, it does not in- 
dicate the direction in which the particle was mov- 
ing (i.e., positive or negative). A technique that 
resolves this directional ambiguity is frequency shift- 
ing; a frequency shift can be chosen to encompass the 
expected range of negative velocities above the zero 
frequency point. 

Fluorometer. Fluorescence intensity is mea- 
sured at the same SV location as the LDV measure- 
ments. (See fig. 1.) The measurement of fluorescent 
dye concentration is based on the fact that the dye 
molecules absorb energy from the incident radiation. 
A molecule that absorbs the input energy is raised 
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to some excited state. While excited molecules dissi- 
pate energy through collisions with other molecules, 
some molecules spontaneously fluoresce as they re- 
turn to their initial ground state. Fluorescence is the 
result of the decay from the excited state to a lower 
molecular energy state. The absorption and emission 
of energy occur at specific wavelengths for each type 
of dye. 

An equation describing the relationship between 
dye concentration and fluorescence (ref. 4) is 

F = $I 0 ( 1 - e~ £bc ) (2) 

where F is the fluorescence intensity, <f> is the quan- 
tum efficiency of the dye, Iq is the incident radiant 
power, e is the molar absorptivity of the dye, b is 
the absorption path length, and c is the molar con- 
centration. The equation shows that an increase in 
concentration increases the fluorescence intensity. At 
very dilute concentration levels, the relationship be- 
comes nearly linear; reference 4 suggests a simplified 
equation 

F = K<f>I 0 £bc (3) 

For a given dye and laser power at a given location 
(no change in b) with K — Constant, equation (3) 
shows a direct relationship between fluorescence in- 
tensity F and concentration c. The linear response 
range of the concentration extends from about 10~ 8 
to 0.1 mg/ml. (See ref. 4.) The dye concentration 
was 0.0245 mg/ml in the injected fluid used to ob- 
tain all of the results reported here. 

The fluorescent dye used in this system was flu- 
orescein. The absorption and emission spectra are 
shown in figure 2. The relative intensity, which is pro- 
portional to dye concentration, is shown as a function 
of light wavelength. The wavelengths from the argon 
ion laser used in the test instrument are indicated at 
the top of the plot. The absorption spectrum rep- 
resents the range of wavelengths which the dye will 
absorb energy from the input radiation. The peak 
absorption occurs at a wavelength of 494 nm and is 
very near the argon laser blue line of 488 nm, which is 
beneficial because nearly maximum energy can be ab- 
sorbed; the emission occurs at the longer wavelength 
of 515 nm. For absorption at any discrete wave- 
length along the absorption spectrum, the dye will 
produce fluorescence throughout the entire emission 
spectrum. If the absorption occurs at a wavelength 
other than at the peak of the absorption spectrum, 
the intensity of fluorescence will decrease across the 
spectrum. 

The peak of the emission spectrum occurs at 
515 nm, which is very close to the argon laser green 


line of 514.5 nm. This produces scattered laser light 
that will contaminate the fluorescence signal; there- 
fore, the fluorescence signal must be filtered both 
electronically and optically. The amplifier for the 
fluorescence signal has a low-pass filter of 30 Hz that 
reduces the contribution of the LDV burst signal in 
the fluorescence signal. However, laser light gen- 
erates broadband shot noise in the photomultiplier 
tube (PMT) some of which will be low frequency 
and will contribute noise to the signal. In addition, 
a 530-nm edge filter that passes light at 530 nm and 
above (50-percent transmission at 530 =t 5 nm and 
less than 1 percent at 515 nm) is used to filter out 
the light below 530 nm that includes the three ar- 
gon laser lines. Therefore, the fluorescence signal 
is derived from light above the edge threshold of 
515 nm. Then, as the concentration of dye in the 
sample volume at the intersection of the two laser 
beams increases or decreases, the fluorescence signal 
level increases or decreases, correspondingly. By in- 
troducing fluorescent dye into one of two interacting 
flows (e.g., a dye-saturated jet emission into a dye- 
free free-stream flow), the mixing of these flows can 
be studied by measuring the variation of fluorescence 
within the flow field. 

The contribution to the fluorescence signal by 
the LDV laser light is optically and electronically 
filtered. However, the fluorescence emission band 
overlaps into the LDV wavelength band. Although 
the fluorescence signal levels are generally lower than 
the LDV levels, an additional broadband shot noise is 
introduced into the LDV signals at the PMT. Narrow 
band-pass filter settings control the contribution of 
the shot noise from the PMT fluorescence signal. 
The low-frequency fluorescence interference caused 
by flow fluctuations was eliminated by 25-MHz high- 
pass frequency filters in the electronic frequency shift 
units (i.e., downmixers) of the LDV system, which 
will be described later. 

Test Facility 

Experiments were performed in the Langley 16- 
by 24-Inch Water Tunnel which is shown in figure 3. 
The tunnel has a useful vertical test section of about 
4.5 ft. The velocity in the test section can be 
varied from 0 to 9 in/sec, which yields unit Reynolds 
numbers from 0 to 7.73 x 10 4 ft -1 , respectively, 
based on a water temperature of 75°F. The typical 
test section velocity is 3 in/sec. Test models in 
the experiments described in this report were sting- 
mounted as shown in figure 4. This system provided 
model rotation in two perpendicular planes with 
ranges of ±15° and ±33°. The center of rotation 
was located on the centerline of the test section 
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and 15.25 in. above the base of the sting. Electric 
motors, which are mounted outside of the test section 
and remotely controlled, drove the model support 
system. Model angles were set with reference to the 
visual indicators within an accuracy of about ±0.25°. 
Other details of the water tunnel are discussed in 
reference 6. 

When the tunnel is full of water, the center of the 
test section is about 10 ft below the water surface. 
The weight of the water causes the test section 
acrylic walls to deflect outward to a maximum of 
approximately 0.090 in. The curvature of the walls, 
which results, alters the alignment of the laser beams 
passing through them. Coincidence of the six beams 
at the same point could not be maintained without 
additional braces; at least two braces were necessary 
to minimize the wall deflection. The braces limited 
optical access and restricted the choices of laser beam 
positions. 

Models 

The first test model was a scaled-down version 
of a wind tunnel propulsion-simulation model with a 
nonaxisymmetric nozzle. Auxiliary water was sup- 
plied to simulate jet flow from the nozzle exit. A 
photograph of this model installed in the test sec- 
tion is shown in figure 5(a); a sketch with princi- 
ple model dimensions is provided in figure 5(b). As 
shown, the axis system origin was centered on the 
nozzle exit. This model was tested at angles of at- 
tack of 0° and 4° and at jet-to-free-stream velocity 
ratios of 0, 1.7, and 3. Dye injection orifices at var- 
ious locations around the aft part of the model pro- 
vided visualization of the flow over the boattail re- 
gion. This dye was introduced into the model from 
the bottom. Fluorescent dye could also be introduced 
into the jet water supply from outside of the tun- 
nel to permit visualization of the jet plume or mea- 
surement of mixing in the jet shear region. The jet 
flow was controlled by a system of valves and a flow- 
meter. Because of the presence of internal dye orifice 
tubes, foam was inserted in the nozzle upstream of 
the throat (fig. 5(b)) to ensure jet flow uniformity. 
The foam tended to compress unevenly and several 
attempts were made to achieve a relatively uniform 
flow distribution. The final test configuration was 
much improved over the model without the foam. A 
drawing of the model with the location of the mea- 
surement planes is shown in figure 6. The planes 
were 15 x 25 arrays (Ay — 0.1 in., Az = 0.0833 in.) 
centered on the nozzle. The measurement plane clos- 
est to the nozzle exit was located one nozzle width 
downstream of the exit (x = 1.032 in.). The second 
plane was located three nozzle widths downstream of 


the exit (x = 3.096 in.). A modified jet model was 
also tested which had a solid constant cross section 
extending one body length downstream from the noz- 
zle exit. Its cross section was the same as the nozzle 
and provided a solid plume simulation of the jet. 

The second test model was an axisymmetric fore- 
body with a strake attached near the nose. This 
model was also a scaled-down representation of a 
wind tunnel model that had been tested to deter- 
mine the effectiveness of strakes for forebody control 
of airplane attitude. The flat strake planform was 
a circular sector located along the model centerline 
with an inboard edge that conformed to the body 
shape. Limited tests were also performed on a cam- 
bered strake that had the same planform shape as the 
flat strake. Dye orifices were located on either side of 
the strake. Figure 7(a) is a photograph that shows 
the model in the test section. The principle model 
dimensions are given in figure 7(b). As shown, the 
camber for the cambered strake was applied longitu- 
dinally; the right edge in the side view is in the X-Y 
plane. Section A- A in figure 7(b) shows the cusped 
cross section of the strake. With the model posi- 
tioned at an angle of attack of 25°, the strake was 
tested at azimuthal angles of 50° and 60° as mea- 
sured from the top of the model. 

Colored dye flow visualization (fig. 7(a)) indicated 
a very coherent vortex shed from the strake; this flow 
was studied with the laser system at three measure- 
ment planes downstream of the strake. The locations 
of the measurement planes relative to the model are 
shown in figure 8. The first plane downstream of the 
nose of the model was located at the downstream end 
of the strake ( x = 3.75 in.) and was an 11 x 25 array 
(Ay — 0.06 in., Az = 0.056 in.). The second and 
third grids (at x = 5.15 in. and x — 10 in.) were 
18 x 25 (Ay — 0.0625 in., Az = 0.0625 in.) and 
17 x 25 (Ay = 0.0625 in., Az = 0.125 in.) arrays, re- 
spectively. All of the measurement grids were located 
with their vertical inboard edges at the model top 
centerline and they extended spanwise in the nega- 
tive Y direction perpendicular to the centerline. The 
vertical extent of the first array was limited by opti- 
cal access between the wall support braces of the test 
section. This also limited the number and location 
of the other measurement grids. 

Experimental Technique 

Optical system setup . The laser measurement 
system operation can be described with the use of 
figure 9. For each measured velocity component, a 
pair of laser beams of the same wavelength cross to 
form a sample volume (SV). The SV is nominally 
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60 /im in diameter and 1.25 mm long and is moved 
throughout the flow field to the measurement loca- 
tion. At each location, three velocity components 
are measured by three pairs of beams at three dif- 
ferent wavelengths. Additionally, fluorescence inten- 
sity is measured from the fluorescent dye, which has 
been introduced into part of the flow field. Mixing 
of the dyed and undyed flows is examined by moving 
the SV through the boundary between the two flows 
(e.g., the shear region between the jet core and free- 
stream flows in fig. 9), and measuring the variations 
in fluorescence intensity. 

The formation of the beam pairs used to measure 
each velocity component is depicted in figure 10. 
Laser light from an argon ion laser is separated 
into its many component colors using prisms and 
mirrors. Three of these colors are selected: violet 
(A = 476.5 nm), blue (A = 488 nm), and green 
(A = 514.5 nm). These beams are split, which forms 
the three pairs of beams. One beam in each pair has 
a Bragg-shift frequency of 40 MHz applied for flow 
direction discrimination. Each beam is sent through 
a fiber-optic cable, passed through a transmitting 
lens outside the test section, and focused at the same 
point to create the sample volume in the test section. 
The scattered light from particulates in the water 
as well as fluorescent light emitted by the dye is 
collected on the opposite side of the test section by 
receiving lenses, which are focused on the SV. 

The collected light is processed as shown in fig- 
ure 11. The data acquisition process begins by focus- 
ing the collected light onto the face of a fiber-optic 
cable. A light-separating optics unit divides the light 
into the three individual laser colors and the fluores- 
cent light color. Each light component is directed 
to a separate photomultiplier tube that converts the 
photons to an electrical signal that contains all of 
the frequency information, which includes the orig- 
inal Bragg-shift frequency of 40 MHz, necessary to 
determine the three velocity components. The elec- 
trical signals are amplified for processing. Each sig- 
nal is high-pass filtered at 25 MHz and its frequency 
downmixed from 40 MHz to 100 kHz, which removes 
the pedestal and reduces the effective Bragg-shift fre- 
quency. This process allows better resolution of the 
Doppler signal frequency relative to the Bragg-shift 
frequency. The signals are then low-pass filtered at 
200 kHz to reduce high-frequency noise before the 
burst counter processing. The velocity resolution of 
the system (smallest velocity measurable) is based 
on the counter resolution, which is about 1.5 percent 
at a velocity of 3 in/sec (^1 mm/sec). The coun- 
ters send frequency information in digital form to the 
data acquisition unit. The fluorescence signal from 


the photomultiplier tube is amplified and sent as an 
analog signal to the data acquisition unit, where it is 
digitized. 

The data acquisition unit collects, multiplexes, 
and buffers the data and then sends it to the com- 
puter. The unit imposes coincidence conditions on 
the LDV data for two or three channels, which en- 
sures that data from the separate channels arrive 
within a user-specified time period. This control is 
necessary for calculating cross correlations such as 
turbulent Reynolds shear stress quantities, for which 
each contributing parameter must be acquired from 
the same measured event (i.e., same particle). In 
addition, the analog fluorescence data is sampled to 
provide data concurrent with the digital LDV data. 
The data from the data acquisition unit is sent via an 
8-bit parallel port to the computer, which is a 32-bit 
desktop system. 

The laser beam transmission and light collection 
system is shown in relation to the water tunnel test 
section in figure 12. The system emits the laser 
beams into the test section from one side; the scat- 
tered light and fluorescent light are collected (off- 
axis) on the opposite side. This arrangement is 
referred to as the forward-scatter mode of light col- 
lection. Two orthogonal beam pairs (violet and blue, 
A = 476.5 nm and 488 nm, respectively) are transmit- 
ted directly on-axis into the test section for measur- 
ing the streamwise and vertical velocities. The third 
beam pair (green, A = 514.5 nm) is transmitted into 
the test section at an angle of 45°, which becomes 
about 32° within the test section after refraction by 
the water. This beam pair allows measurement of the 
velocity component normal to the beam angle bisec- 
tor, which is then used to calculate the on-axis lateral 
velocity component in the test section coordinate sys- 
tem. Figure 13 is a photograph of the transmitting 
and receiving optics. 

The transmitting and receiving optics are po- 
sitioned by a traverse control system that can be 
either operated independently of the computer in 
manual mode or controlled interactively by the com- 
puter system. The positioning system controls seven 
motor-encoder units: three X-Y-Z units on the 

transmitting side, an auxiliary unit for the off- 
axis laser beams, and three X-Y-Z units on the 
receiving side. The encoder resolution for the aux- 
iliary motor-encoder system is 0.000025 in/pulse 
(0.635 jzm/pulse). The resolution for the other 
six motor-encoder units is 0.00005 in/pulse 
(1.27 (i m/pulse). 

The characteristics of the sample volume for 
each of the three component colors are given in 
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Table I. LDV Sample Volume Characteristics 


[ See figure 14 ] 



W/d, in./in., of — 

Characteristics 

0/18.2 

4/15.2 

8/12.2 

12/9.2 

16/6.2 


Axial: A = 476.5 nm (violet), dj = 3.85 pm 


/, in 

19.4 

20.4 

21.4 

22.4 : 

23.4 

dm , pm ... . 

51 

54 

56 

59 

62 

Ijji , mm .... 

1.10 

1.16 

1.21 : 

1.27 

1.32 


Vertical: A = 488 nm (blue), dj — 3.94 p m 


/, in 

19.4 

20.4 

21.4 

22.4 

23.4 

dm : pm .... 

52 

55 

58 

60 1 

63 

l m mm .... 

1.12 

1.18 

1.24 

1.30 

1.36 


Lateral: A = 514.5 nm (green), dj — 4.98 pm 


/, in 

19.8 

22.1 

24.4 

26.6 

28.9 

dm: pm .... 

56 

63 

69 

76 

82 

Z m , mm .... 

1.45 

1.62 

1.78 

1.95 

2.12 


table 1 based on the distances symbolically defined in 
figure 14. Note that when the LDV system is run in 
the coincidence mode (i.e., all three velocity measure- 
ments are taken from the same particle), the effective 
measurement region is the intersection of the sample 
volumes formed by the various beams. Because the 
sample volumes are approximated by ellipsoids (i.e., 
the locus of the beams where the power is e“ 2 of 
the peak of the Gaussian beam power distribution), 
the intersection of the SV’s can be approximated by 
calculating the length of the intersection of the max- 
imum ellipsoid dimensions. For the beams in the 
center of the test section, this calculation yields an 
intersection length of about 220 g m. This represents 
the effective SV length for measurements made in 
the coincidence mode because it is within the focus 
length of the receiving optics. The maximum inter- 
section diameter is about 56 gm, which is equal to 
the smallest single SV diameter. 

Laser light sheet . Flow visualization was ob- 
tained by laser light sheet excitation of the fluorescein 
dye. The laser light sheet was generated by sweeping 
an argon ion laser beam of 488 nm at a 500-Hz rate. 
This gave a light sheet with an essentially uniform 
power distribution. Details on the laser light sheet 
system can be found in reference 7. The laser light 
sheet was oriented perpendicular to the longitudinal 
body axis. Flow patterns visualized by the fluores- 
cence were photographed end-on from downstream of 
the model by reflecting the image to the outside of 


the tunnel with a mirror placed at an angle of 45° in 
the bottom of the test section. 

Seeding material . The seeding particles used 
in the experiments can be described with figure 15. 
The top of the figure shows a scanning electron 
microscope photograph of the seeding particles in a 
water sample taken from the water tunnel during a 
test run; the bottom of the figure shows a particle size 
distribution plot. The number of particles increases 
as the average particle size decreases; the smallest 
particle size range shown is 5 to 10 gm. The average 
size for the sample was 13 gm. The analysis indicated 
that only 0.5 percent of the particles were greater 
than 100 gm. The particles were identified as being 
mostly hydrated iron oxide, Fe 203 (i.e., rust), with 
an average density of 3.6 gm/cm 3 . By analytical 
methods of references 8 and 9, the majority of the 
test seeding particles were determined to be small 
enough to follow the flow accurately. 

Model test conditions and setup . All tests 
were run at a free-stream velocity of 3 in/sec. The 
nonaxisymmetric propulsion model was tested at an- 
gles of attack of 0° and 4°. The jet-to- free- stream 
velocity ratios were 0, 1.7, and 3. In addition to the 
planar data, vertical velocity profiles were measured 
near the nozzle exit and axial scans were taken down- 
stream of the exit on the model centerline. 

The axisymmetric forebody with a strake was 
tested at an angle of attack of 25°. The baseline 
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flat strake was tested at strake azimuthal angles of 
( ft = 50° and 60° relative to the top of the model. 
The cambered strake was also tested at the same two 
angles. However, because of limited optical access to 
this model, only one F-axis scan for each azimuthal 
angle was made. 

At the beginning of each test, the LDV system 
optics were aligned in an empty test section. The 
model was then inserted in the tunnel test section 
and the tunnel was filled with water. After the water 
was circulated to minimize any density gradients 
because of temperature variations, the tunnel flow 
was set to the test velocity of 3 in/sec. The axis 
system origin for these data was centered on the 
nozzle exit for the jet test and at the nose of the 
forebody for the strake test. 

Measurements were initiated; as the beams were 
moved from point to point, the burst counter data 
rates were monitored to ensure that optical alignment 
was maintained. At each point, up to 1000 samples 
were acquired with a maximum acquisition time of 
1 to 2 min. All data reported herein were calculated 
from the samples taken at each point. 

Analysis of Results 

The results of the experiments are presented 
in this section for the nonaxisymmetric propul- 
sion model (jet test) followed by the axisymmetric 
forebody-strake model (strake test). An estimate of 
the uncertainties in the calculated results are given 
in tables AI(a)-AI(c) of the appendix. 

Nonaxisymmetric Propulsion Model 

Flow visualization . Laser light sheet flow vi- 
sualization results are shown in figures 16 and 17. 
The qualitative turbulent development of the jet at 
a jet-to-free-stream velocity ratio of 1.7 and angle of 
attack of 4° is shown in the sequence of photographs 
in figure 16. The nozzle flow was visualized with the 
laser light sheet perpendicular to the jet longitudi- 
nal axis at the exit (x/W = 0), one nozzle width 
downstream of the exit ( x/W — 1), and three noz- 
zle widths downstream of the exit (x/W = 3). Dye 
traces in the jet flow show the increased mixing at 
the periphery of the jet as the laser light sheet was 
moved downstream. The flow from the boat tail re- 
gion is shown in the sequence of photographs in fig- 
ure 17 for the same flow conditions as figure 16 but at 
different downstream locations of x/W. In figure 17, 
only the dye external to the model is shown, which 
mixed with the jet in a very short distance from the 
exit. 


Mean velocity results . The axial velocity u 
is the component of primary interest in jet stud- 
ies. Average axial velocity data taken for the jet at 
several spanwise stations are shown in figure 18 at 
a = 0° and a jet-to-free-stream velocity ratio of 1.7. 
The uncertainty in this quantity was estimated to be 
about =1=2.4 percent of [Too* The maximum velocity 
occurred at the station y = 0.3 in. and decreased 
on either side. (The offset in the +F direction was 
most likely caused by some residual nonuniformity 
of the foam in the nozzle throat.) In addition, the 
velocity profiles at the closest measurement station 
x — 1.032 in., (fig. 18(a)) showed a velocity deficit 
between z — 0 and —0.2 in. for two scan locations 
(y — 0 in. and —0.3 in.). This was probably caused 
by either the model support strut and the associ- 
ated wake momentum deficit or affected by the flow- 
conditioning foam in the throat. Comparison of the 
two axial measurement stations reveals a reduction 
in peak velocities at the aft station x = 3.096 in. 
(See fig. 18(b).) The velocity reduction effect of in- 
creasing longitudinal distance is shown directly for 
each of the jet-to-free-stream velocity ratios in fig- 
ure 19. In addition, the effect of velocity ratio at 
each station is shown by comparing the velocity lev- 
els in figures 19(a) and 19(b). 

Flow characteristics that are not discernible from 
line plots are often revealed by examining data con- 
tours in the measurement plane. Contours of axial 
velocity at a = 0° are shown in figure 20(a) at the 
first measurement station downstream of the nozzle 
exit (x/W = 1). The axial velocities were maxi- 
mum in the irregularly shaped region near the center 
of the measurement grid. A ring of low-momentum 
flow around this high-speed region that resulted from 
flow separation on the boattail is shown. Also, a re- 
gion of reduced axial velocity is evident in the lower 
center section of the grid. This velocity defect was a 
result of the momentum deficit in the wake shed by 
the model support strut. The effect of angle of at- 
tack is shown in figure 20(b). The effect was evident 
in two concentrated low-velocity regions in the top 
corners of the grid. They represent the reduced axial 
velocity regions associated with body-shed vortices 
from the model upper surface corners. 

Vertical distribution of axial velocity . A 
limited number of surveys were made near the nozzle 
exit on the leeward side of the model. The surveys 
were made outside of the jet in a vertical direction at 
angles of attack of 0° and 4°. Vertical distributions 
of axial velocity at the nozzle exit are shown in 
figure 21 for three jet-to-free-stream velocity ratios. 
(Note that the bottom of the vertical axis is not at 
zero; however, the vertical upper extent of the nozzle 
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is shown in fig. 21. The vertical upper extent of 
the body is 0.615 in.) An increase in jet velocity 
ratio increased local axial velocity u for a substantial 
range of Z. This is indicative of jet interference 
effects on the nozzle boattail. A significant region of 
reversed flow existed with the jet off. A comparison 
of velocity profiles which includes a solid-plume jet 
simulator is shown in figure 22. The solid-plume 
jet simulator was a constant-dimension extension 
aft of the nozzle exit. The model with the solid 
plume generated more reverse flow than the jet-off 
configuration, which indicated a trend opposite to 
what would be expected. The solid-plume jet model 
was expected to represent the flow better than the 
jet-off configuration. This result indicates that the 
jet flow must be provided for proper simulation of 
off-body flow physics. A direct comparison of the 
effect of angle of attack at three velocity ratios is 
shown in figure 23. In all three cases, angle of attack 
increased the local axial velocity u profiles. 

Axial velocity decay . The jets in this study 
are effectively three-dimensional incompressible tur- 
bulent jets. In general, incompressible turbulent jets 
are characterized by the presence of three distinct re- 
gions in axial velocity decay along the jet axis. (See 
refs. 10 and 11.) These regions are illustrated in fig- 
ure 24 and are classified as follows: 

1. Potential core (PC) region . This region is the 
section from the nozzle exit at point 0 to point A. 
Mixing that occurs at the jet boundaries has not 
yet permeated the entire flow field and leaves a 
region that is characterized by a uniform velocity 
on the axis equal to the jet exit velocity. 

2. Characteristic decay (CD) region. This region is 
the section from point A to point B. Within this 
region, the axial velocity decay is dependent upon 
nozzle exit geometry; velocity profiles in the plane 
of the minor axis of the exit are found to be similar 
whereas, those in the plane of the major axis 
are nonsimilar. Therefore, this region is termed 
characteristic of the initial geometry. 

3. Axisymmetric decay (AD) region. This region is 
the section from point B to point C where point C 
is not a specific, well-defined location. The axial 
velocity u decays in this region in accordance with 
u oc x -1 , which is characteristic of axisymmetric 
jets. The entire flow approaches axisymmetry 
and thus exit geometry is no longer a factor. 
Velocity profiles in this region are similar in both 
planes of symmetry. Far downstream beyond 
point C, a fourth region of fully axisymmetric flow 
is observed. 


Again, note that most of the characteristics of the 
jet in the PC and CD regions (i.e., near field) appear 
to be determined completely by the exit geometry 
and that flow in the AD and fully axisymmetric 
region (i.e., far field) is independent of the exit 
geometry. 

In accordance with the discussion by Rajaratnam 
(ref. 11), the test model produces a bluff compound 
jet. The exit aspect ratio of the model in this 
study (width/height of W/2bo = 1.9) qualifies this 
jet to be bluff as well as compound because of the 
motion of the surrounding fluid. Bluff jets originating 
from square nozzles have been studied by Trentacoste 
and Sforza. (See ref, 12.) When the characteristic 
decay region is excluded, they discovered that the 
potential core of a square jet is followed very quickly 
by an axisymmetric decay region. These trends were 
apparent in the current experimental results. 

Experimental observations of the current config- 
uration of bluff compound jets were made at jet- 
to-free-stream velocity ratios of 1.7 and 3.0, with a 
nozzle exit aspect ratio ^1.9. The axial scans were 
made in the near field up to x/2b$ ~ 9 (far field is 
considered to be x/2bo > 100), where x is the axial 
distance and &o is the nozzle half-height. Figure 25 
presents the axial velocity decay of the jet for both 
jet velocity ratios. Despite the scatter in the data for 
Uj/Uoo = 3, the potential core (u/u maLX — 1) seems 
to extend to x/bo — 2 for both jet velocity ratios. The 
uncertainty in u/u m ax was ±0.9 and ±0.5 percent for 
Uj/Uoo — 1-7 and 3, respectively. The similar behav- 
ior for both jet velocity ratios supports the hypoth- 
esis made by Sforza (ref. 10) that the characteristics 
of the PC region are determined by the nozzle exit 
geometry and not the initial velocity. Following the 
potential core region, the two curves decay linearly 
on the log- log plot. (The curve for the velocity ratio 
of 3 follows the trend better than for 1.7.) The jet 
of this investigation seems to follow the trend of a 
bluff jet, which is characterized by the potential core 
followed by an axisymmetric decay region. 

Turbulence results . The turbulence results 

were calculated on the basis of the statistical analysis 
technique of mean and standard deviation for experi- 
mental samples. Note that other sources of apparent 
turbulence such as vortex meandering, which could 
contribute to deviations from the mean flow, were 
not addressed here. 

The Reynolds stresses reported herein will actu- 
ally be the stresses per unit density. Two exampl es of 
actual Rey nolds normal and shear stresses are pu f u ' 
and pw f v! , respectively. Reynolds normal stresses 
are analogous to pressures because they represent the 
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reaction of a surface to turbulent flow of momentum 
(e.g., pu f2 dS) through the surface in the direction of 
the momentum. The distri bution of average axial 
Reynolds normal stress u'u' is shown in figure 26 at 
two measur emen t stations. The uncertainty in this 
normalized u f u f was ±0.067 percent. The levels of 
this stress decreased with increased distance from the 
nozzle exit and with reduction of velocity ratio. The 
asymmetry in the profile is accentuated at the higher 
velocity ratio. In general, the depression in the pro- 
file at the station x = 1.032 in. for the lower velocity 
ratio is characteristic of jets because the velocity fluc- 
tuations tend to be significant in the shear layer at 
the periphery of the jet. The vertical distribut ions 
of the average transverse Reynolds normal stress v f v f 
and the vertical Reynolds stress w f w f are shown in 
figures 27 and 28, respectively. The estimated un- 
certainties for these quantities were ±1.5 percent and 
±0. 042 percent, respectively. Trends similar to those 
for v!u ' occur with the change in veloc ity ratio and 
distance. Howeve r, th e le vels in v f v f were signifi- 
cantly higher than u f u f and w f w f levels. Other exper- 
imental studies (e.g., Hinze, ref. 13) show variations 
from complete isotropy in the se st resses to domina- 
tion by the axial compo nent v!v! over the smaller 
vertical vj 'w' and lateral v f v f components. The high 
indicated v f v f levels in this study probably result pri- 
marily from a larger background transverse turbu- 
lence intensity (i.e., rms velocity fluctuation) vf in 
the water tunnel test section and, to a lesser extent, 
from the propagation of er rors through the coordi- 
nate transformation of vfvf . 

After the model tests, the velocity fluctuation v f 
was measured in the empty test section. The off- 
axis optical equipment used to measure the trans- 
verse velocity was dismantled and remounted to send 
the green laser beams of 514.5 nm through a perpen- 
dicular tunnel sidewall to measure the transverse ve- 
locity component directly. The resulting orthogonal 
system revealed a nonisotropic turbulence environ- 
ment. The directly measured transverse component 
v f was consistently 15 to 25 percent greater than vf , 
which was consistently 30 to 44 percent greater than 
it/. Possibly one contribution to the large v f is the 
transverse oscillation generated upstream where the 
water encounters the back wall and turns the cor- 
ner without the aid of turning vanes. Another pos- 
sible contribution is the capture of the tunnel-free 
surface waves which are drawn down into the test 
section by the curved lip at the upper back wall. 
Finally, without an isolation system, acoustic dis- 
turbances from the pump or from cavitation on the 
back side of the speed-control butterfly valve could 
be transmitted through the water into the test sec- 


tion in some predisposed direction. These oscillations 
could also cause localized tunnel wall movement that 
would alter the instrument focus and produce appar- 
ent turbulence. 

The distribut ions of a verage Reynolds shear 
stresses u f v f , u f w ( , and w ( u f are shown in figures 29 
31, respectively, for two measurement planes and two 
velocity ratios. These tangential stresses are anal- 
ogous to shear stresses because they cause a shear 
stress reaction on the surface. They represent cross 
correlations of turbulent fluctuating velocities and 
show the turbulent transport of momentum in an or- 
thogonal direction. For example, the u f v f stresses in 
figure 29 represent the turbulent transport of X mo- 
mentum in the Y direction. In general, the u f v f corre- 
lations were positive; therefore, a positive increase in 
v! corresponded to a positive increase in v f . The es- 
tima ted u ncertainty in this stress was ±0.22 percent. 
The v f w ; stresses in figure 30 represent the turbulent 
transport of Y momentum in the Z direction. The es- 
timated uncertainty in this stress w as ±0 .10 percent. 
These stresses were lower than the u f v f stresses but 
showed a strong reversal near the centerline in both 
measurement planes, especially for the high-velocity 
ratio. Note that the value of the shear stress, whether 
positive or negative, indicates the relative am ount 
of momentum transfer due to turbulence. The w'u f 
stresses in figure 31 represent the transport of Z 
momentum in the X direction. T he e stimated un- 
certainty was ±0.036 percent. The w f u' stresses were 
similar in magnitude to the u l v f stresses; however, 
the character of the distribution was similar to the 
vfw 1 stresses. The distributions changed sign on ei- 
ther side of the jet centerline which is characteristic 
for these shear stresses in a jet. (Note that the ef- 
fective center of the jet is not at z = 0.) A contour 
plot of the w f u f stress in figure 32 shows the three- 
dimensional nature of the distribution in the jet. 

Fluorescence results . The distribution of the 
local average, nondimensionalized fluorescent dye 
concentration across the jet is given for two mea- 
surement locations in figure 33. The estimated level 
of uncertainty was ±0.71 percent. Dye was injected 
into the jet flow and the fluorescence intensity mea- 
surement showed the dispersion and extent of mixing 
with the free-stream flow. As shown in the figure, the 
spread of the jet increased with increasing distance 
downstream of the exit. The dimensional peak con- 
centration value c was reduced at the aft measure- 
ment station, which indicated that the spread of the 
dye as a result of mixing reduced the concentration 
at the center. 


10 


The fluctuations (i.e., standard deviations) in flu- 
orescent dye concentration are shown in figure 34. 
The estimate of uncertainty was ±0.57 percent. 
These fluctuations occur where fluid that is tagged 
with dye mixes with undyed free-stream fluid. High 
fluctuation values indicate high levels of mixing as 
occurs in a shear layer. For a jet, the fluctuations 
would be high in the shear layer surrounding the jet 
and be near zero in the free stream (i.e,, no dye) 
and in the center of the jet where the dye concentra- 
tion remains constant (i.e., no mixing). The distri- 
butions in figure 34 indicate the increased mixing in 
the downstream direction as the shear layer grows. 
The centerline value increased as the turbulent mix- 
ing reached the center of the jet. Contours of the 
concentration fluctuations are shown in figure 35 for 
the first measurement plane downstream of the noz- 
zle exit. The contours map out the shear layer be- 
tween the free stream and the jet by identifying the 
region of mixing. 

The cross correlation between velocity and con- 
centration fluctuations is indicated by calculation of 
the statistical covariance of these quantities similar 
to the way in which the mean shear stresses were 
calculated. For the streamwise component u f 

N 

£ (Ui - u)(Ci - c ) 


where Ui and c\ are instantaneous values. Equa- 
tion (4) is the numerator (i.e., the covariance) of 
a correlation coefficient for u ! and d. Therefore, 
if this quantity is zero, v! and d are uncorrelated 
and independent. These relationships are shown in 
figures 36-38. All of the distributions as well as the 
concentration fluctuations exhibited a reduction at 
the center of the jet. The correlations between u f 
and d are indicated in figure 36 and were similar 
for the two velocity ratios at each measurement sta- 
tion. The estimated uncertainty, when normalized 
by Uqo and the maximum c value for the scan, was 
±0.015 percent. The values were positive-correlated 
in the jet shear regions and indicated a stronger re- 
lationship between the axial and concentration fluc- 
tuations than did the other two velocity fluctuations. 
This implies that the axial fluctuations were the pre- 
dominant mechanism in shear layer mixing. The v f d 
relationship is indicated in figure 37 and the w f d 
cross correlations are shown i n fig ure 38; their lev- 
els are clearly lower than the u ( d quantities. Their 
respective estimates of uncertainty were ±0.056 and 
±0.005 percent. 


Axisymmetric Forebody With Strake 

Flow visualization results on the axisymmetric 
forebody indicated that the strake position of 50° 
provided the most coherent vortex structure. A 
position of 60° was also investigated to provide a 
comparison with the strake configuration of 50°. 
(Only these two strake angles were investigated in 
this study.) Observations confirmed the confinement 
of fluorescent dye in the vortex core both visually 
and quantitatively in the fluorescence measurements, 
which negated the possibility of turbulent mixing in 
the vortex. The highly three-dimensional nature of 
the vortex provides a flow for which the LDV is well 
suited and one that is very different from the jet flow 
shown previously. In the flow visualization figures, 
the flow is moving from right to left in each figure 
and the nose of the body is in the left side of each 
figure. 

Flow visualization . Flow visualization results 
are shown in figures 39 and 40 for strake azimuth 
angles of 0 — 50° (baseline) and 60°, respectively. 
The flow visualization technique used was laser light 
sheet excitation of fluorescein dye with visualization 
from the downstream perspective. The development 
of the strake vortex in the axial downstream direc- 
tion was visualized by positioning the laser light sheet 
perpendicular to the body longitudinal axis. The 
four stations shown represent the midstrake position 
(fig. 39(a)), followed by the three measurement sta- 
tions where quantitive data were taken. The develop- 
ment of the strake vortex for 4> = 50° indicated that 
the vortex remained laminar along the entire length 
of the model. Figures 39(b)-39(d) show the develop- 
ment of a body vortex shed from the model surface 
inboard of the strake vortex. This is characteristic 
of slender bodies at incidence. Both vortices were 
rotating in the counterclockwise direction in the fig- 
ures. Initially, the vortex generated by the strake at 
60° (fig. 40(a)) showed a similar characteristic to the 
vortex generated by the strake at 50°. However, as 
the vortex migrated downstream, the appearance of 
the streaklines from the strake at 60° was less well- 
defined and the vortex appeared to be larger. This 
would imply a vortex was formed with a broader ve- 
locity distribution and a larger axial velocity deficit; 
examination of the mean velocity distributions shows 
this to be true. 

Mean velocity results . Average u, u, and w 
velocity distributions across the vortex of the strake 
at 50° are given in figures 41(a)-41(c), respectively, 
at several vertical z stations. The respective un- 
certainties, when normalized by U 0 c , were estimated 
to be 0.011 or -0.012, 0.031 or -0.033, and 0.0085 
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or —0.0098. The scans were taken at the first mea- 
surement station x = 3.75 in. downstream of the 
body nose. The location of the vortex center was 
z = 1.35 in. and y — 0.725 in.; the peak axial u 
(fig. 41(a)) and tangential (vertical, w in fig. 41(c)) 
velocities occurred at this Z location. The distribu- 
tions of transverse velocity v (fig. 41(b)) indicated 
that the center of the vortex (z = 1.35 in.) was 
the location of minimum transverse velocity. This 
is reasonable because there should be minimal ve- 
locity through the center of the vortex along the 
scan T-axis. The axial velocity distributions in fig- 
ure 41(a) revealed that velocity levels were different 
on either side of the vortex as a result of the presence 
of the flow field around the body. The centerline of 
the body was at y = 0; z increased away from the 
body in the vertical direction. The axial velocity 
was greater on the body side of the vortex, perhaps, 
because of the accelerated flow around the side of 
the body. One exception was the closest 2 station 
(z = 1.08 in.) where the flow had apparently sepa- 
rated and the measurements were made in the wake 
of the body. The axial velocities from the other scan 
locations show a gradual decrease from y = —0.4 
to 0, although still larger than Uoo • This may result 
from the gradual approach to the wake region where 
the effect of the characteristic momentum deficit is 
increasing. Also, there was a curious bump in the 
distributions outboard (y « -1.0) of the main vortex 
axial velocity deficit. This distortion is the result of 
the presence of the strake wake (i.e., separated free- 
shear layer) before it wrapped into the strake vortex. 
A similar bump was present in the tangential (verti- 
cal) velocity distributions in figure 41(c). 

A comparison of the three mean velocity distri- 
butions at the first measurement station for the two 
strake roll angles is shown in figures 42(a)-42(c). The 
axial velocity distributions are plotted in figure 42(a). 
The distribution for the strake at 60° is displaced in 
y because of the physical displacement of the tip of 
the strake when rotated outboard an extra 10°. The 
axial velocity distributions reveal that the strake at 
60° generated a larger axial velocity deficit than the 
strake at 50°. In addition, the width of the velocity 
deficit is slightly broader. The transverse (fig. 42(b)) 
velocities v do not show significant differences in 
terms of levels. The vertical (fig. 42(c)) velocities 
w show slightly higher peak values and a greater dis- 
tance between upper and lower velocity peaks for the 
strake at 60°. This result indicates that the vortex 
from the strake at 60° was spread out more as was 
implied by the flow visualization shown earlier. In 
general, as vortices approach breakdown, the ratio of 
peak tangential to axial velocities increases. In this 


investigation, therefore, the vortex from the strake 
at 60° would appear to reach breakdown first. 

The effect of cambering the strake at 50° is shown 
for the three mean velocities in figures 43(a)-43(c). 
The X location for the cambered strake was slightly 
different because of the limited optical access which 
resulted from the braces around the test section. 
The cambered strake showed much less axial veloc- 
ity deficit than the baseline flat strake in figure 43(a), 
which indicated that the baseline strake vortex would 
likely break down sooner. No significant difference 
in transverse velocity was noted in figure 43(b). The 
vertical velocities in figure 43(c) indicated that the 
cambered strake generated slightly higher velocity 
peaks. The width of the vortex viscous subcores, 
defined as the region between the peaks of the tan- 
gential velocity distribution, appeared to be quite 
similar. 

Contours of mean axial velocity for the strakes at 
50° (baseline) and 60° are shown in figures 44 and 45, 
respectively. Comparison of the minimum values in 
the center of the contours shows the greater extent of 
the velocity deficit for the strake at 60°, which again 
indicates the tendency of the vortices from the strake 
at 60° to break down sooner. 

Streamwise vorticity and circulation . 

Streamwise vorticity was calculated from the aver- 
age v and w velocity data. Identification of regions of 
concentrated vortical flow enables location of vortex 
centers and the determination of relative local vortex 
strengths. Contours of streamwise vorticity in the 
three measurement planes are shown for the strakes 
at 50° and 60° in figures 46 and 47, respectively. The 
estimated maximum uncertainty in streamwise vor- 
ticity was calculated to be 2.59 sec -1 . This rather 
large value is based on the worst-case vertical veloc- 
ity gradient, which exists in the center of the vortex. 
In other regions of the vortex flow, the uncertainty 
would be greatly reduced. The peak vorticity gen- 
erated by the strake at 50° was stronger at the first 
measurement station than for the strake at 60°. How- 
ever, the vorticity from the strake at 50° decreased 
more at the second station than did the vorticity from 
the strake at 60°, which was then slightly higher. 
The peak vorticity at the final measurement station 
was about the same for both strakes. Evidence of 
the wake from the strake before it is wrapped into 
the vortex exists on the lower right-hand side of the 
measurement grids. 

The variation of streamwise vorticity across the 
center of the vortex for the three measurement sta- 
tions is shown in figure 48 for the strake at 50°. 
The vorticity peaks were reduced at the downstream 
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measurement stations. This reduction was accompa- 
nied by a slight lateral spread in the vorticity distri- 
bution on the right-hand side of the peaks. However, 
the vorticity appears to be decreasing overall in the 
downstream direction probably because of viscous 
dissipation, especially, when considering the flow low 
test Reynolds number of 2.5 x 10 4 ft -1 . The effect of 
the strake wake was observed as a bump in the previ- 
ously mentioned velocity distributions. There is also 
a bump from the effect of the strake wake in the vor- 
ticity distributions on the — Y side of the vorticity 
peak. 

The amount of vorticity contained in the viscous 
subcore has been examined. The streamwise vor- 
ticity calculated at the grid points is compared to 
the vertical velocity distribution in figure 49 for the 
strake at 50° at the second measurement station. 
On the right-hand side of the figure, a very small 
amount of vorticity appears inboard of the point 
(y « “0.56 in.) where the peak negative velocity 
begins to move back toward zero. On the outboard 
side of the figure (y « —1.0 in.), the same trend was 
apparent. However, the wake of the strake added a 
small bump of vorticity before it diminished. Over- 
all, most of the vorticity appeared to exist in the area 
defined as the vortex subcore. The remainder of the 
vorticity is contained in the so-called rotational core 
(i.e., outer core) which is very small for this trail- 
ing vortex when compared to the outer cores of wing 
leading edge vortices. (See ref. 14.) This outer core 
behavior is consistent with the theory that the area 
outside of the outer core can be represented as a po- 
tential vortex which implies irrotational (i.e., non- 
vortical) flow. 

Calculation of total circulation indicates total vor- 
tex strength and allows vortex comparisons. Circula- 
tion values were calculated for the two strake angles 
at each of the measurement stations. Because of the 
sign convention, circulation on the side of the model 
where the measurements were taken was calculated 
as a negative value. At the first station, the vorticity 
distributions were cut off by limited optical access, 
especially for the strake at 60°, as shown in figures 46 
and 47. This prevents any comparison at the first 
station. At the second station, circulation values for 
the st rakes at 50° and 60° were —2.76 in 2 /sec and 
—3.03 in 2 /sec, respectively, with an estimated un- 
certainty of 0.012 in 2 /sec. Integration over the en- 
tire grid at the third station included some additional 
negative vorticity from other sources, especially for 
the strake at 50°. By integrating over that part of the 
grid which excluded this additional vorticity, circula- 
tion values of —2.94 in 2 /sec and —3.13 in 2 /sec were 
calculated for the strakes at 50° and 60°, respectively. 


Comparison of these circulation values implies that 
the vortex from the strake at 60° was stronger and 
therefore likely to breakdown sooner. 

Fluorescence results . The spanwise distribu- 
tion of normalized fluorescent dye concentration is 
shown in figure 50 for the strake at 50° along the 
vortex center scan line at the three measurement sta- 
tions. The estimated uncertainty in this quantity was 
calculated to be ±0.93 percent. Baseline values have 
been removed from c and c max - The dye remained 
confined to the viscous core and no turbulent mixing 
was evident. 

Conclusions 

An instrument that simultaneously measures 
three velocity components and mixing (from flu- 
orescence intensity) in a flow field by application 
of laser Doppler velocimetry and fluorometry was 
demonstrated in the Langley 16- by 24-Inch Water 
Tunnel. Comprehensive three-dimensional flow-field 
data were acquired from two different types of flow. 
First, a jet flow generated by a propulsion simulation 
model with a nonaxisymmetric nozzle was investi- 
gated. Flow-field data were measured at two stations 
downstream of the model exit. Mean axial velocity 
profiles were presented to show the general charac- 
teristics of the jet at two jet-to-free-stream velocity 
ratios and two angles of attack. Vertical profiles of 
axial velocity above a solid-plume jet simulator re- 
vealed more reverse flow than for the jet-off configu- 
ration. This indicates that jet flow must be provided 
to properly simulate off- body flow physics. 

Six turbulent Reynolds stresses were calculated 
from the measured data. The Reynolds normal 
stresses were maximum near the center of the jet; 
the Reynolds shear stresses were generally greater in 
the jet shear regions. Dye concentration fluctuations 
were mapped in the shear region which indicated 
the region of free-stream and jet flow mixing. The 
covariance of velocity fluctuations with fluorescent 
dye concentration fluctuations indicated that the 
predominant velocity fluctuation in shear-layer mix- 
ing was the longitudinal velocity fluctuation. How- 
ever, all three velocity fluctuations contributed to 
shear-layer mixing. From examination of the three- 
dimensional characteristics of the flow, the three- 
dimensional velocity-field data and the simultaneous 
predominant mechanisms in mixing were determined. 
This had not been accomplished previously. 

The experimental results from a test of an axisym- 
metric forebody with a strake revealed the character- 
istics of a laminar vortex flow with an axial velocity 
deficit. The flow field generated by a strake at 50° 
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was compared with that at 60°. The results indi- 
cated less axial velocity deficit for the strake at 50°. 
Streamwise vorticity and circulation were calculated 
from the measured data. Comparisons between the 
two strake angles indicated that the vortex generated 
by the strake at 60° was stronger and may break 
down sooner than the vortex from the strake at 50°. 
These results are important for the determination of 
optimal forebody-strake configurations. 

The quality of the measurement system has been 
studied by performing an uncertainty analysis of the 
experimental results. The calculated uncertainties 
were different for the results of the two tests. In 
general, uncertainties in the mean velocities varied 
between 1 and 7 percent of free-stream velocity. 


Fluctuating velocities varied between 1 and 6 percent 
of reference values. The uncertainties in mean lateral 
velocity v were consistently greater than for the other 
two mean velocities. The uncertainties in the lateral 
velocity fluctuations ?/ were also generally greater 
than for the other two components. An independent 
direct measurement of v f indicated a consistently 
greater value than for either v! or it/ in the free 
stream. Also, the consistently greater v and v f 
uncertainties resulted from the effect of resolving the 
measured off-axis components into the orthogonal 
coordinate system. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
July 13, 1994 
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Appendix 

Analysis of Experimental Uncertainties 

The estimation of uncertainties in the results from measurements made by the instrument was based 
on standard principles outlined in the ANSI/ASME Measurement Uncertainty. Part I-~Instruments and 
Apparatus . (See ref. 15.) A current book by Coleman and Steele (ref. 16) that explains the contents of 
the document was used in the present analysis. The concepts of bias and precision errors are defined as fixed 
and random errors, respectively, that occur in experimentation. The bias errors are the fixed, systematic, 
or constant errors that induce an offset from the true value of the quantity being measured. The precision 
errors are random variation or repeatability errors of the quantity being measured. The total uncertainty is 
the root-sum-square (RSS) combination of the bias and precision error estimates and is given by 

u RSS =(£ 2 + ^ 2 ) 1/2 

where B is the bias error estimate and P x is the precision limit. The precision limit is the product of the precision 
error estimate (precision index) and £95 from the statistical t distribution that gives a 95 -percent confidence 
level to the estimate of the random precision error. For experimental sample sizes greater than 30, reference 15 
recommends that £95 — 2 for 95-percent confidence estimates. The equation for the total uncertainty becomes 

^RSS = [S 2 + (2S) 2 ] 1/2 (Al) 


where S is the precision error estimate or precision index. 

Coleman and Steele (ref. 16) describe the propagation of uncertainty in a general uncertainty analysis. If a 
result r is a function of many variables xi, # 2 , # 3 ? • • x^y, the uncertainty U r in a calculated result is 


U r = 



The Ux s are the uncertainties in the measured variables X{. In a detailed uncertainty analysis, the bias and 
precision limit estimates must be propagated separately. Then the estimates are combined into the total 
uncertainty as in the RSS method in this case. If the bias and precision limits in the measurements of the 
different variables are not independent of each other, cross terms exist in the expressions for the bias and 
precision limits, respectively, as follows: 


B r — 
P r = 



dr dr 
dxi dxj 


B Xi B x . p B '' ( 1 — 6ij ) 


dr dr 
dxi dxj 




(A2) 


where the pij ’s are correlation coefficients associated with each type of error and 6ij is one if i = j and zero if 
i 7 ^ j. If the error terms are independent, the p^’s are zero and the cross terms dissapear from equations (A 2 ). 

The bias limit for r is 


B r — 



and the precision limit for r is 


where the Bx s and the P x ’s are the bias and precision limits for the measured variables X{. The tabulated results 
show the propagated precision error estimates (i.e., precision indices) instead of the precision limit. However, 
the total uncertainty contains the precision limit P — t$§S. The equation for the propagated precision error 
estimate is 


S r 




+ ••■ + 


dr 

dxjsj 


?X N 



where the S x ’ s are the estimated precision errors in the measured results. 


The data acquisition program for the instrument calculates the untransformed results from the measured 
quantities. These results are modified by a coordinate transformation to yield final results. Therefore, the 
uncertainty analysis was required to take into account this transformation. The propagation of the uncertainties 
into the calculated results, which included coordinate transformation, was performed by a method similar to 
that used by Morrison et al. (See ref. 17.) 

Categories of errors for the estimation of untransformed result uncertainties are given in the next section. 
The development of coordinate transformation equations are then described. Finally, the method for calculating 
uncertainties in the final transformed results will be shown. All calculated quantities were normalized by 
appropriate terms as shown in the tables at the end of this appendix. Instead of point by point, overall test 
total uncertainties were estimated for each calculated quantity. Worst-case conditions were generally used, 
which provided an estimated upper limit on uncertainties. Final estimated uncertainties in tables AI(a)-AI(c) 
are presented with two significant figure accuracy because additional significant figures would be inappropriate 
for such estimates. 


Uncertainties in Results Calculated From Measured Quantities 

Four quantities are measured by the instrument: Ai/ W jV , A v u ^ g , A v u ^ and which are the three frequency 

components (yields three velocity components) and the fluorescence intensity in the laser Doppler velocimeter 
(LDV) coordinate system, respectively. From these untransformed measurements, the following 17 quantities 

are calculated: u , F, FJ; a/?/ 2 , V^ 2 , \Ai/ 2 ; u f u f , v f v r , w f w f \ u f v f , v f w f , w f u f \ \ 


v f B, 


and 


c, c . 


The elemental uncertainties in the quantities measured by the velocimeter were estimated by the methods of 
Patrick (ref. 9) and Meyers. (See refs. 18-20.) Three categories of errors were determined for the velocimeter: 
laser beam geometrical errors, processor errors, and seeding-induced errors. These errors are itemized in 
table AII(a). Elemental processor errors were estimated for the fluorometer. Finally, statistical errors (i.e., 
precision errors) were calculated and are listed in tables AII(a)-AII(c) for the untransformed results. The 
following sections address the various categories of error sources. 


Velocimeter Errors 

Laser beam geometrical errors . Seven sources of error were considered in this category. They included 
positioning uncertainty, beam orientation, crossbeam angle, finite probe volume, negative velocity, incomplete 
signal, and frequency- broadening biases. The last three errors were either not significant or not applicable to 
the system. See reference 9 for a more detailed explanation. 

Processor errors . Six sources of error were examined in this category. These included clock synchroniza- 
tion, quantization, electronic noise, comparator tolerance, threshold limit, and pedestal removal filter. The last 
three sources were insignificant. 

Seed-induced errors . Four types of errors were considered in this category. They included seed particle 
lag, velocity bias, flow distortion, and Bragg bias. Only the first two were significant in this study and velocity 
bias was only significant in the strake test. No corrections were made for velocity bias. 


Fluorometer Errors 

The principal fluorometer errors are attributed to the analog-to-digital (A/D) processor which digitizes 
the fluorescence signals and the amplifier. The A/D converter which changes the amplified analog fluorescence 
signal to a digital signal is a 16-bit converter with a ±5-V range. According to reference 16, the accuracy, which 
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implies bias, is ±V 2 of the least significant bit (LSB). This is equivalent to one half of the digital resolution 
times the voltage range or 

AjF=± K^) =±0 ' 0763mV 

The maximum fluorescence signal processor output voltage F max was about 80 mV. Therefore, the error was 
estimated as A F/F max « ±0.1 percent. The maximum fluctuations in dye concentrations were about 35 percent 
of F max (~0.028 V) so that the error estimated in the fluctuations was ±0.27 percent (0.0763 mV/0.028 V). 
The accuracy specified by the amplifier manufacturer was ±0.1 percent of the signal. The total errors estimated 
for the fluorescence quantities were calculated by the RSS of the individual A/D converter and amplifier error 
estimates. The estimates are summarized in table AII(c). 

Statistical Errors 

Statistical errors occur because of the averaging nature of data processing in the LDV system. Finite 
numbers of samples are averaged at each measurement location. As a result, the true standard deviation, 
which indicates the level of precision error for an infinite number of samples, cannot be known. The estimate 
of the standard deviation is defined as 


jfTT - s > 2 

1=1 


1/2 


where 


1 N 
u ^ E 

i = 1 


and where N is the number of samples, U{ is the ith velocity sample, and u is the mean velocity at the 
measurement location. 


If several sets of mean velocities u or mean dye concentrations c are calculated, the sample distribution of the 
mean values is assumed to follow a normal distribution about the parent population mean value. Therefore, the 
standard deviation of the normally distributed means represents the precision error in the mean (i.e., standard 
error of the mean) and was estimated as 



or 


S c — 


Vn 


Because the mean velocity samples are assumed to be normally distributed, the standard deviations (and 
thus approximately the rms velocities) are assumed statistically distributed in a chi-square manner. (See 
ref. 9.) However, for large sample sizes (N > 50), the chi-square distribution approaches a normal distribution. 
Therefore, the precision error in the rms turbulence velocities can be estimated as the standard error in the 
standard deviation (refs. 9 and 21) by 


S u ' ~ 



V2N 


or 


S<f = 


V2N 


Patrick (ref. 9) also assumed that the Reynolds shear stresses were chi-square distributed and for large 
sample sizes (JV > 50), the precision error could be estimated as 


S- 7-7 = 


u'v' 

V2N 
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The rms velocity-concentration fluctuation correlations were treated similarly to the Reynolds shear stresses 
and the precision errors were estimated as 

Zp 

V ~ V2N 

In this study, a maximum of 1000 samples were taken. At times, data rates were lower and time did not 
permit the acquisition of 1000 samples. Therefore, 800 samples were used in the precision error estimation 
calculations. 

All of these estimates are for quantities in the LDV coordinate system. However, because the LDV system 
has a third off-axis laser beam component, the estimates were transformed into the water tunnel orthogonal 
coordinate system. This was accomplished by propagating the results with the coordinate transformation 
matrices, which are derived in the next section. 

Uncertainties in Final Calculated Results 

The component of velocity u g is measured with an off-axis set of laser beams. A coordinate transformation 
was used to calculate the final results in the tunnel coordinates. Thus, any quantity calculated with the 
u g component will be affected by this transformation— including uncertainty estimates. Therefore, the 
transformation equations will be derived in general and then applied to this particular system. 


Coordinate Transformation 

The derivation of the transformation equations follows the development of Morrison et al. (See ref. 17.) 
In figure Al, the axes shown as dark lines form an orthogonal X-Y-Z coordinate system. The total velocity 
vector defined in this system would be 

U = u\ -f vj -f- w\a 


where i, j, and k represent unit vectors in the A, Y , and Z directions, respectively. This represents the 
velocity defined in the tunnel coordinate system. The velocity defined in the LDV coordinate system would be 
represented generally as 

Uldv — u v e v 4- u g e g + 


where e v ,e g , and are the unit vectors in the LDV coordinate system directions. In general, the two coordinate 
systems have the same origin but the axes are not coincident. The relationship between the respective unit 
vectors can be represented by 



<* 12 <*13 

<*22 <*23 

<*32 <*33 



and, therefore, the relationship between the velocity vectors is 


u v | 

Oil 

012 

013 

r « 

u 9 f = 

021 

022 

023 

\ v 

u b ) 

.031 

032 

033. 

( w 


Because the measurements were made in the LDV 
defined by 


u 

v 

w 


cii 

C 21 

C31 


coordinate system, the inverse relationship was used as 


02 03 

02 03 

<%2 03 



(A3) 


where 

foil = KP 1 


The derivation of the aij matrix will be followed by the derivation of the inverse c\j matrix. 

One example of an LDV coordinate axis is shown by the light line in figure A 1(a). Here, 6 g is the angle 
in the X-Y plane between the X-axis and the e g unit vector in the LDV coordinate system. The other 
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two orthogonal LDV components are shown in figure Al(b), rotated about the K-axis by an angle 6 V in the 
X-Z plane. The angles 9 V and 9 5 are the angles between the water tunnel X-axis and the measured e v and 
unit vectors, respectively; the angle 9 ^ will always be 9 V 4- 90° from the X-axis. Based on the angles between 
the tunnel X-axis and the LDV-measured velocity vectors, equations were derived for the corresponding unit 
vectors. An example of the derivation geometry is shown in figure Al(c) from which the unit vector was 
derived. The equations are 

e v = i cos 9 V - f k sin 9 V 


The resulting equation for 


e g = i cos 9 g — j sin 9 g 

e\y = i cos 9fy 4- k sin 9 ^ 
velocities in terms of the a^j matrix is 


U V ] 

I cos 9 V 

0 

sin 9 V 

u 9 

> = COS 9 g 

I 

5' 

os 

'O 

0 

u-b J 

! cos 9b 

0 

sin 9b 


From the fact that the product of a matrix and its inverse equals the identity matrix, the terms in the c t j 
matrix can be determined by the equation 


cos 9 V 

0 

sin 9 V 

\c n 

c 12 

c 13 

cos 9g 

— sin 6 g 

0 

C 21 

c 22 

c 23 

cos 9\y 

0 

sin 9ft 

L C 31 

C32 

C33 



'1 

0 

0 ‘ 

— 

0 

1 

0 


0 

0 

1 


Multiplication of the terms to form nine equations results in the q ? - solution matrix 


sin#/, 

sin (0 b -6 v ) 


- sin 9 V 
sin (6 b -6 v ) 


sin $ b cos 0 g —i — sin 9 V cos 0 g 

smdgSmiOb-Ov) sin & g s'm(0^-0 v ) sin 0 g 

— cos 9\y q cos 6y 

sin {9 b -9 v ) sin {9 h -~9 v ) 

The above inverse matrix C{j was verified by multiplying it by the original a ij matrix to yield the identity 
matrix. 

The transformation equations for the correlations between the fluctuating velocities u' and the fluorescence 
fluctuations c f were similar to those for the mean velocities. They are represented in condensed form as 


\ v c __ [cij] u g c 
ly^'J [tz'c' 

From the definition of the Reynolds normal and shear stresses, 


N 

£ UiUi 


u'u f = 


N 

E UiVi 


/ f i~ 1 

u'v f = — — uv 


where the UiS and are total instantaneous velocity values of the sum of the mean velocity u or v and the 
corresponding instantaneous fluctuating velocity u f or v f . Substitution of the matrix equations (A3) for the 
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Ui and u values into equation (A4) yields the following transformation equation for the u f u f Reynolds normal 
stress: 


N ( c ll c ll u v,i u v,i H - c ll c 12 u v,i u g,i T c ll c l 3 u v,i u b,i 

+C12CH Ug,iU v ,i + C\2Ci2UgiUg^i + Ci2Ci3UgjU bi 
i=l +C!3Ciiu b jU v j + ci3Ci2U b jUgj + cizCisUb'jUbj) 

N 


{c\\Cl\U v U v + C\lC\ 2 UvUg + C\\C\^U v Ub 
+ Cl2Cl\UgU V + C\2C\2UgUg + Ci2C\^U g Ub 

+C\ZC\\UbU v + Ci3Ci2U b Ug 4- C13 C !3*W>) 


Similar equations were derived for the other two normal stresses and the three shear stresses but are not 
shown here. The end result was a set of transformation equations that is expressed as the following matrix 
equation: 


u f u f 


'CllCll 

c 12 c 12 

C13C13 

2cnCi2 

2cnc]3 

2C12C13 


u v u v 

v f v f 


c 21 c 21 

c 22 c 22 

c 23 c 23 

2C21C22 

2C21C23 

2C22C23 


u 'g u 'g 

w f w f 


c 31 c 31 

C32C32 

C33C33 

2C31C32 

2C31C33 

2C32C33 


u 'b U 'b 

u'v' 


C11C21 

C12C22 

c 13 c 23 

C11C22 + C12C21 

C11C23 + C13O21 

c 12 c 23 + C13C22 


u' v u'g 

v f w f 


C21C31 

C22C32 

^23 c 33 

C21C32 + C22C31 

C21C33 + C23C31 

c 22 c 33 + C23C32 


u' v u' b 

_ w f u f _ 


-C31C11 

C32C12 

c 33 c 13 

c 31 c 12 + C32CH 

^31 c 13 + C33C11 

C32C13 + C33C12- 


l u 'g u bi 


For the LDV coordinate system under study, the following quantities are defined: 


e v = o° 

e b = 90° 

Q g = 32.037° (inH 2 0) 


This leads to the definitions of the two transformation matrices as 


[ c ij\ ~ 


cos 0 ( 


1 __ 


and of the Reynolds normal and shear stresses as 


1 

0 

0 

0 

0 

0 • 

( cos 0,J \ 2 

( -1 f 

0 

—2 COS Og 

0 

0 

y sin 0g ) 

V sin Og J 

sin 2 0 g 

0 

0 

1 

0 

0 

0 

COS Og 

0 

0 

-1 

0 

0 

sin 0,i 

sin&g 

0 

0 

0 

0 

COS 0g 

-1 

sin8g 

sin 9g 

0 

0 

0 

0 

1 

0 . 
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The equations for the transformed results are 


U — Uy 


V = u v 


COS 0 g 

sin 6 g sin 6 g 


W = Uf } 


(A6) 


.1 n> 


COS Og U'gC 

qy = V C — — — — 

v sin 6 0 sin 9 n 


w'd = u' b d 


(A7) 


1 TT 2 cos °9 ~TT 


, 7 / COS 0g\2 — — — 

” ” “ ( sin 6 g ) UvUv + sin 2 6 g ® sin 2 


it/u/ = 


(A8) 


where 


/cos0p 

vfv r — — — -r~ \uLvL, — 


sin e 9 J~ v v sin 6 g v g 


cos 6 t 


- 'u'ui — 


si nOgJ^ tinOgW 


W f U f — U f b U f b 


^rms 


= \fu r 


(A9) 


^rms = V v'v' 


^rms 




Of course, because the fluorescent dye concentration c and fluctuation d are scalar quantities, they are 
transformed without modification to the final results. 

Application of Uncertainty Propagation Into Calculated Results 

By reference to the transformation equations (A6) for the mean velocities, the equations for the propagation 
of the bias limits and precision indices (i.e., precision errors) can be derived. Because 


U — Uy 

du 


8Uy 


= 1 


w = qi} } 
dw 


dub 


= 1 


and, for example, 


By = 


\2 

dUy^ 


1/2 


21 


then the propagated bias limits and precision errors for the u and w mean velocities are equal to those for the 
measured quantities. That is, 


B u — B Uv S u — S Uv B w — B Ub S w — S Ub 
For the on-axis velocity component v, the bias equation from equations (A2) becomes 
2 / dv „ \2 ( dv „ \2 ( dv n Y dv dv n _ 

~ ( duT Buv ) + ( P>~ Bu S ) + ( fift, Bf) <> ) + rfTTT. ffnZ Pu ^ U 9 Bu v Bu g H 






< 9 tln 


The bias in in this study was ±0.1° (±0.0017 rad). The errors in the measurements of u v , u g , and 9 g are 
assumed to be independent. The cross terms drop out (p Uv u g = Pu g 9 g ~ Pu v 6 g = 0), which results in 


By = 


cos 

sin 


Br, 


+ 


1 \2 / 9 \ 9 o 11/2 

" sin <9 + \ Uv CSC ° 9 + Ug CSC ^ COt ® 9 ) (■ A1 °) 


Similarly, for the precision error, 


Sy = 


SSf^ 1 +| - 


sin 6 a B Ug 


There is no precision error in 9 g ; the error in 9 g is a bias error. 
As before, the total uncertainty is then calculated as 


Uu = 


Uv = 


UnJ = 


Bl + (2S w f 
Bl + (2 S„) 2 
5|+(25^) 2 


1/2 

1/2 

1/2 


11/2 


The bias limits and precision errors for the Reynolds normal stresses were calculated from the transformation 
equations (A8). Again, because 


'W ! W f = u b U b 


the bias limits and precision estimates were derived as 


B— = B-rr S— = S~ B— = B— 

u'u IhiUt, U r u r Ui,U-n w W U b U b 


IXytly W UyUy 

From the formulas for the Reynolds stresses, 


5— = 5— 

urw' u b u b 


B j , — B , , — B i f — 0 ) 


XLylLy 


HgUg 


bb 


B—j—r — B—tt ~ B—j—y = 0 

'U'v'U'g ^Q^b 


(All) 


This implies that the errors in the Reynolds stresses are strictly statistical precision errors before they are 
transformed into the tunnel coordinate system. 

The errors in the measurements of u f v v! v and u l v u f g or u g u f g and u f v u g may not be independent which would 
generate cross terms in the error equations. However, the extent of dependence is not known; also, the 
magnitude and sign of the correlation coefficients of the cross terms are not known. Therefore, the measurement 
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errors in u f v u f v , u g u g , u f v u g , and are assumed to be independent for this quantity, which eliminates cross 
terms. The bias limit in v f v f can be estimated from 


B- 


^Lb-^ Y + ( Lb — Y + 


V / dv f v f Y 

du f v u g u ' vU 9 J y dd g & 9 J 


From the transformation equations (A8) for the v f v f Reynolds normal stress, the bias limit is 


B\j = 

V V 


cos6 A 2 b 

sin On ) u v u v 


+ 


B-j ~r + 1 


sin 2 6 g u g u 9 


—2 cos 0 g 
sin 2 0 o 


B-t-t 


+ 


-2 cos 0 g 
sin 3 On 


(u ( v u f v + Ugii'^j + 2 u f v u f g esc 0 g ( 2 cot 2 0 g + 1) 


Bi 


From equations (All), the first three terms are zero and the bias is a function only of the error in 0 g . The 
precision index has a similar form of 


v v 


COS Og 
sin 0 n 


5-7- 


+ 1 


1 

sin 2 0, 


] + 


-2 cos 0 g 

sin 2 On u v u g j 


where S-i—r is the measured shear stress precision error estimate. 

v! v v! g 

The transformed fluctuating (rms) velocities are calculated as the square root of the transformed Reynolds 
normal stresses. Therefore, 


u 1 ^ — u Tms 


= V Vtv! = \j UyUy — 


U Vn 


and 


— w Tms 


- \J u 'b u 'b - 


u b„ 


Similar to the mean velocities, the bias and precision errors are 


Bu rms Bu Vrms B w rms ^ u b Tms ^rms ^t'rms ^rms ^ u b rms 

transformation equations (A8) and by noting that u 2 rms = u f v u f v , and ^ 2 rms = u g u g > the fluctuating 

ms is 

/= f / ens \2 „ 1 „ 2cos0„—— 1 1/2 


Urms = vVl/ = ( U 2 V + 

\ sin 0 g ) rms sin" 1 

the bias limit is derived from 


2 cos Og 

sin 2 6 g U9rms 


• 2 a 
sm z 0 g 


r 2 _ / dvrms R \ 2 
" rms \du Vims 


' dv rl 


, f - ^rms 


’rms ry V 
% Be >) 


^rms r> V , f ^rms R \ 

^ + \dv\y g «) \ d6 g 

appropriate terms in the above equation and algebraic manipulation yield 
7 COS 0 fJ y U Vrms E I 2 ( Bu grms U 9, mS ¥ . ( ~ COS 8g V 

A sin Og) Wrms “ t,rms J ysin 2 9 g v Tms J ^ sin 2 0 g v ims J 

no / ^ \ o 


B l 

‘'rms 


+ ( U Ls + “Ls) + U >9 csc M 2 COt2 °9 + l ) 

,erm is zero. 
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The precision errors are 


2 _ \( COsflg Y 


U V n 


The uncertainty analysis for the Reynolds shear stresses was based on the transformation equations (A9). In 
this case, the velocity m easu reme nts fr om the same particle were required for correlation of velocity fluctuations. 
Therefore, each of the u f v f and v'w f stresses had a third cross term in their expressions for bias limits and 
precision errors. The errors were assumed to be perfectly positive-correlated and the correlation coefficients 
were set equal to one. The equations for the bias limits are as follows: 


v! v v! v esc 2 9 g + v! v v! g esc 6 g cot 6 g j 2 Bq 


+ 





ll v Ug 


= (~ U 'v U b csc2 6 9 + u 9 u b CSC e 9 cot e 9f B \ 


(COSO^B \2 

\ sin 0 g u' v u' b ) 


+ 


-B —\ 2 

U 9 U b \ 

sin On 


/ cos^ Y - 1 V B 

V sin 0„ J\ sin 6„ ) u W b 


Vg y ailic/p J \ Olll Ug / \OH lVg, 

By application of equations (All), only the first term in each equation is nonzero. Also, 


B zrzr - B zrj = 0 


a v u b 


The precision errors are calculated as follows: 


u'v y sin 0 g u vu v ) 


'-STrrrV 


+ 


uv lL g 


sin On 


( cos^y^LY 5 

V sin On J\ sin 0 g J u W g 


-S-r-r\2 


2 /cos 0, V . ( cosflg Y zMc ? 

v'W V sin On uWj \ sin On I + ^sin 0 g A sin 0 g J W g u' b 


and 


*5 — S , , 


Fluorescent dye concentration c is not measured directly but is calculated from the following equation for 
the fluorescence intensity given previously: 

F = K<f>I 0 ebc 

where b represents the size of the sample volume and the other terms are as defined before. The solution for 
the concentration is 

F 

C " K$I 0 £b 

The bias limit for the concentration is calculated from 


B r = 


( dc \2 ( dc „ v (dc „ v 

\dF BF ) + \di 0 Bl °) + \db B 7 


1/2 
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with the assumption that the measurements of F, Jo, and b are independent and the contributions from K, <J>, 
and e can be neglected. The resulting equation for the bias limit is 


B c = c 



211/2 


For the tests discussed in this paper, the variations in laser intensity were assumed to be small enough to be 
ignored during the sample period. Also, the sample volume was assumed to be completely filled with dye during 
the sample period and to not vary in size. This results in an equivalent relationship between the fractional 
bias limits for concentration and fluorescence intensity. From the above assumptions, the relative bias in the 
concentration results is proportional to the relative bias in the fluorescence measurements and either 


or 


B c 


B c = 


Bp 

F 

Bp 


K<f>I 0 £b 

Because the correlations of velocity and dye concentration fluctuations were transformed exactly as the 
mean velocities were, their transformation equations are similar to equations (A7). From the formulas for the 
correlations, the bias e rrors in t he u f d and w f d correlations are zero. With the assumption of independent 
measurement errors in u f v d , u g d and 6 g , the resulting equations for bias limits are 


B—rj = B 

u'c 


v! v d 


- B w’d - B u'd ~ 0 


and 



( cos p V , 

\ sin Og <d) + 



sin 6 g 


+ 


u f v d esc 2 9 g + u g d esc 6 g cot 



where only the last term is nonzero. Similarly, for the precision errors, 






and 





An example calculation of the positive bias error in v for the jet test is presented here, 
equation (A10) for this error is 


Recall that 


By = 


+ + (■“” csc2 e ’ + csc 


All velocities used in this equation were nondimensionalized by U 0 G . The angle 6 g = 32.037° with a bias of 
±0.1° (±0.0017 rad). Values of u v /Uoo — 1.5 and u g /Uoo = 0.05 were used. The bias in u^/U 0 0 is 0.0207 
and —0.0213 from table AII(a) and includes contributions primarily from laser beam geometry and particle 
lag errors; the bias in u^/U 0 0 is 0.0216 and —0.0221 from the same sources of error. The positive bias errors 
in Imj/Uqo and u^/Uoo result in the three squared terms of equation (A10) being equal to 0.00109, 0.00166, 
and 0.00008, respectively. By RSS methods, the total positive bias error in v/Uoo is 0.0532. Clearly, the first 
two terms contribute the most to this error because they are predominately influenced by the error in the 
crossbeam angle measurement. This result is typical of most systems and greater accuracy in crossbeam angle 
measurement will have the greatest effect in reducing the bias error. 

This concludes the discussion of the propagation of uncertainty for the basic calculated quantities. All of 
the transformed final results are summarized in tables AI(a)-AI(c). All transformed quantities were normalized 
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by the appropriate terms as shown in the tables. The next section describes the use of the total uncertainties 
in mean velocities in the calculation of the uncertainties in streamwise vorticity and circulation from the strake 
test. 


Estimates of Uncertainty in Vorticity and Circulation 

The calculations of streamwise vorticity and circulation were performed with the results of the study of the 
axisymmetric forebody with a strake. The derivations of equations for vorticity and circulation were based 
on references 22-24. The following general equations are used to formulate uncertainty estimates for the 
calculations of streamwise vorticity and circulation. 


Errors in vorticity calculations. The general equation for the streamwise vorticity and the discretized 
approximation is 


/ dw 
V ~dy 



A w Av 
Ay A z 


The total estimated uncertainties in velocities were the combined bias and precision error. They were 
combined into an alternate bias error for calculation of the error in vorticity. With the assumption that 
v and w , and y and 0 were independently measured, an analysis based on total uncertainty values of the 
components of vorticity yields the error 


du x p 

” \ J-> A w 


2 

+ 



B Av 

~Kz 


2 


The following equation was derived for the error in streamwise vorticity: 


B 







To obtain the worst-case estimates of uncertainty, the minimum spacing of Ay and Az for the grid was used. 
The largest absolute values for the uncertainties in the v and w mean velocities were used. The calculations 
yielded the following result for the uncertainty in streamwise vorticity: 




OJ x 


6.687 + 0.00028 
2.59 sec -1 


Clearly, the first-order terms in Ay and Az are the dominant terms and the effect of the error in measurement 
position is small for this calculation. Because the uncertainties in mean velocities are greatest in high-velocity 
gradient regions, the level of error in u) x is probably less in lower uj x regions where velocity gradients diminish. 
(See fig. 49.) 


Error in circulation calculations. Circulation was calculated by two methods: a line integral of the 
velocity around a closed path (i.e., rectangular cell) and an area integral of the vorticity, which resulted in 
circulation values that were within 0.04 percent. The equations and their approximations for each method, 
respectively, are 


r 

r 


<p v dy + w dz « ^ Vi Ay + Wi Az 
C * =1 

J ui x dS « uj ce \\ Ay Az 
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The error in the calculated circulation for the line integral method can be estimated as follows: 


2=1 

Bl = 4[(A yB v f + ( vB Ay ) 2 + (A z B w ) 2 + ( wB Az ) 2 ] 

From typical peak values of v and w and the values defined in the previous section, the following uncertainty 
in the circulation was calculated: 

By = 0.0124 in 2 / sec 

The error in the calculated circulation for the area integral method can be estimated as follows: 




2 (dv „ v / dv „ v ( a r _ v 

-Bp = (Ay A,z B^) 2 + (u x Az B Ay ) 2 + (u x Ay B Az ) 2 

From a peak vorticity value and other values previously defined, the uncertainty in circulation was estimated 
as 

By = 0.0088 in 2 / sec 


Analysis of Uncertainty Results 

The significance of the final uncertainty results can be examined by comparing them with reference values. 
The reference value for all fluctuating quantities was 0.53f/oo and was determined by calculating the averages of 
the mean plus standard deviations of the velocity fluctuations u\ 1 /, and w f from the jet test. The uncertainties 
in the mean values are already referenced to the free-stream velocity and are shown in table AI(a). The 
remaining uncertainty results relative to 0.53C/oo (or (0.53f/oc) 2 for stress terms) are given as percentages in 
table AIII. The results from the strake test are significantly lower because of much lower turbulence values 
(ssll percent of Uqo)- The velocity-concentration quantities were referenced to 0.53t/oo and 0.1 c ma x (c m ax = 
45 mV). 

All quantities associated with the on-axis v component (including mean values in table AI(a)) had the 
greatest uncertainties because of the resolution of the off-axis geometry into the on-axis ^-component direction. 
This system has been analyzed by Neti and Clark, Yanta and Ausherman, and Bell, Rodman, and Mehta. (See 
refs. 25, 26, and 27, respectively.) They formulated a relation for the error in the on-axis velocity component 
as a function of the off-axis angle, an assumed fractional uncertainty in measured components, and the ratio 
of on-axis velocity to measured velocity. When the velocity ratio is small, the errors in on-axis velocity can 
be extremely large (>100 percent). This implies that the errors estimated on a point-by-point basis may be 
considerably greater or less than the values estimated in this study. However, the error values given in the 
tables are a reasonable estimate of the orders of magnitude. 
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Table AI. Uncertainty Estimates for Transformed Results 
(a) Mean and fluctuating (rms) velocities, Uqq = 3 in/sec 


Error 

U 

Uoc 

V 

W 

\fvfi 

u x 

\Zj/2 

Uoc 

\/ u/2 

Jet test 

Bias 

+0.0207 

-0.0213 

+0.0532 

-0.0545 

+0.0129 

-0.0138 

±0.00504 

±0.0154 

±0.00389 

Precision 

±0.00580 

±0.0198 

±0.00456 

±0.00410 

±0.0150 

±0.00322 

Total 

uncertainty 

±0.024 

+0.066 

-0.067 

±0.016 

±0.0096 

±0.034 

±0.0075 

Strake test 

Bias 

+0.0101 

-0.0112 

+0.0254 

-0.0281 

+0.00822 

-0.00957 

±0.00195 

±0.0119 

±0.00441 

Precision 

±0.00168 

±0.00888 

±0.00110 

±0.00119 

±0.00600 

±0.000775 

Total 

uncertainty 

+0.011 

-0.012 

+0.031 

-0.033 

+0.0085 

-0.0098 

±0.0031 

±0.017 

±0.0047 


(b) Reynolds normal stresses, Uqc = 3 in/sec 


Error 

vJv! 

Ux 

v'v' 

U% 

w , w > 

W 

Jet test 

Bias 

0 

±0.00027 

0 

Precision 

±0.00067 

±0.00728 

±0.00042 

Total 

uncertainty 

±0.00067 

±0.015 

±0.00042 

Strake test 

Bias 

0 

±0.00005 

0 

Precision 

±0.00006 

±0.00149 

±0.00002 

Total 

uncertainty 

±0.00006 

±0.0030 

1 

±0.00002 


(c) Reynolds shear stresses and fluorescence- related data, U 0 0 = 3 in/sec 


Error 

vfv f 

TJT 

u oo 


w f u ' 

u'd 

v’d 

w'd 

c 

c' 


W 

t^oc c max 

PocCmax 

Poo c max 

c max 

c max 

Jet test 

Bias 

±0.0000728 

±0.0000622 

0 

0 

±0.0000238 

0 

±0.0014 

±0.0029 

Precision 

±0.00110 

±0.00050 

±0.00036 

±0.00015 

±0.00028 

±0.00005 

±0.00346 

±0.00244 

Total 

uncertainty 

±0.0022 

±0.0010 

±0.00036 

±0.00015 

±0.00056 

±0.00005 

±0.0071 

±0.0057 

Strake test 

Bias 

±0.0000029 

±0.0000012 1 

0 




±0.0014 

±0.0029 

Precision 

±0.00012 

±0.000066 

±0.00002 




±0.00461 

±0.00326 

Total 

uncertainty 

±0.00024 

±0.00013 

±0.00002 




±0.0093 

±0.0071 
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Table AIL Uncertainty Estimates for Untransformed Results 
(a) Mean and fluctuating (rms) velocities, U 0 c = 3 in/sec 

[Jet test] 


Error 


ft 

ft 

U a c 

~uzr 

\J u 'b 

Beam geometry 

Position 

uncertainty 

±0.00364 


^0 

±0.00020 

±0.00031 

±0.00015 

Beam 

orientation 

±0.00197 

±0.0128 

±0.0128 

±0.00154 

±0.00116 

±0.00168 

Crossbeam 

angle 

±0.0195 

±0.0174 

±0.00130 

±0.00107 

±0.00170 

±0.00084 

Finite probe 
volume 

±0.00667 

0 

0 

±0.00467 

±0.00708 

±0.00338 

Processor errors 

Clock 

synchronization 

0.0000075 

0.0000019 

0.0000025 

0.0000037 

0.0000013 

0.0000010 

Quantization 

^0 

«0 



«0 

«0 

Electronic 

noise 

0.000155 

0.000262 

0.000062 

0.0000016 

0.0000057 

0.0000002 

Seed-induced errors 

Particle lag 

-0.0049 

-0.0049 

-0.0049 

0 

0 

0 

Velocity bias 1 








Totals 

Total bias 

±0.0207 

-0.0213 

±0.0216 

-0.0221 

±0.0129 

-0.0132 

±0.00504 

±0.00738 

±0.00389 

Precision 

errors 

±0.00580 

±0.00927 

±0.00456 

±0.00410 

±0.00655 

±0.00322 

Total 

uncertainty 

±0.0237 

-0.0242 

±0.0285 

-0.0288 

±0.0158 

-0.0165 

±0.00962 

±0.0150 

±0.00752 


^ot computed. 
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Table AIL Continued 


(a) Concluded 
[ Strake test ] 


Error 

Uy 

Uoo 

Un 

£ 

T7ST 

\f u 'g 2 

A 2 

Beam geometry 

Position 

uncertainty 


±0.000133 

±0.000177 

±0.000067 

±0.000240 

±0.000160 

Beam 

orientation 

±0.00513 

±0.00513 

±0.00394 

±0.000623 

±0.000310 

±0.00069 

Crossbeam 

angle 

±0.00833 

±0.00785 

±0.00654 

±0.000309 

±0.00083 

±0.000202 

Finite probe 
volume 

±0.0025 

±0.00362 

±0.00304 

±0.00182 

±0.00652 

±0.00435 

Processor errors 

Clock 

synchronization 

0.0000025 

0.0000019 

0.0000025 

0.0000002 

0.0000004 

0.0000001 

Quantization 

«0 


«0 

«0 

«0 

«0 

Electronic 

noise 

0.000155 

0.000262 

0.000062 

0.0000002 

0.0000016 

«0 

Seed-induced errors 

Particle lag 

-0.0049 

-0.0049 

-0.0049 

0 

0 

0 

Velocity bias 2 







Totals 

Total bias 

±0.0101 

-0.0112 

±0.0101 

-0.0112 

±0.00822 

-0.00957 

±0.00195 

±0.00658 

±0.00441 

Precision 

errors 

±0.00168 

±0.00449 

±0.00110 

±0.00119 

±0.00318 

±0.000755 

Total 

uncertainty 

±0.0106 

-0.0117 

±0.0136 

-0.0144 

±0.00851 

-0.00982 

±0.00308 

±0.00915 

±0.00467 


2 Not computed. 
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Table AIL Concluded 


(b) Reynolds normal stresses, U 0 0 = 3 in/sec 


Error 

uW v 

Uoo 

u ’q u ’g 

P 

u oc 

Jet test 

Bias 

0 

0 

0 

Precision 

±0.00067 

±0.00171 

±0.00042 

Total 

uncertainty 

±0.00067 

±0.00171 

±0.00042 

Strake test 

Bias 

0 

0 

0 

Precision 

±0.00006 

±0.00040 

±0.00002 

Total 

uncertainty 

±0.00006 

±0.00040 

±0.00002 


(c) Reynolds shear stresses and fluorescence-related data, Uoo = 3 in/sec 


Error 

UyUg 

7/2 

u oc 

u 'i u 'b 

W 

t'ST 

u'yCf 

u' g d 

v! h d 

c 

c' 

t'oo c max 

t'ooQnax 

Poo c max 

c max 

c max 

Jet test 

Bias 

0 

0 

0 

0 

HDHH 

0 


Ell— 

Precision 

±0.00060 

±0.00016 

±0.00036 

±0.00015 





Total 

uncertainty 

±0.00060 

±0.00016 

±0.00036 



±0.00005 

±0.00706 

±0.00568 

Strake test 

Bias 

0 

0 

0 




±0.0014 

±0.0029 

Precision 

±0.00007 

±0.00004 

±0.00002 




±0.00461 

±0.00326 

Total 

uncertainty 

±0.00007 

±0.00004 

±0.00002 




±0.00933 

±0.00714 
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Table AIII. Final Estimated Total Uncertainties in Fluctuating 
Quantities Relative to Reference Values 



c / c max 


1 


1 
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u n — Av djr 
In which: 

u n Velocity component of particle 
normal to laser beam bisector 

Av Doppler-shift frequency of scattered light 

dj Fringe spacing 

Figure 1. Dual beam LDV and fluorescence probes. 
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Forebody coordinates, in. 
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Cambered strake 

3-percent biconvex wrapped around 
NACA-63 mean line with 
planar outer edge 


Top view Side view 



(enlarged) 


(b) Principle dimensions. All dimensions are in inches. 
Figure 7. Concluded. 









Collection of 
scattered light 



Computer 


Figure 11. Data acquisition system. 

















Water tunnel 

Traverse and ontics, sec ^ on Traverse and optics. 



Figure 12. Water tunnel laser instrument. 
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Figure 13. Laser instrument installed at Langley 16- by 24-Inch Water Tunnel. 
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Indices of refraction 
n i = 1 .0 (air) 
n 2 = 1.43 (Lucite™)* 
«3 = 1.333 (water) 


Lucite: trademark of E. I. du Pont de Nemours & Co. 
Figure 14. Beam refraction through media. 





(a) Scanning electron microscope photograph of seeding. 


Thousands 
per liter 

300 -i 



Average size, pm 


(b) Particle size distribution. 

Figure 15. Water tunnel seeding for laser Doppler velocimeter (natural particulates). 



(a) Scanning electron microscope photograph of seeding 


Thousands 


Average size, pm 
(b) Particle size distribution. 

Figure 15. Water tunnel seeding for laser Doppler velocimeter (natural particulates) 








(b) x/W = 1. 



= 1.7. 
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jc, in. 

1.032 

3.096 



Figure 19. Effect of measurement location downstream from nozzle exit on axial velocity profile of non 
axisymmetric boattail model, y — 0 in.; a = 0°. 





I x = 1.032 in.; Uj/Uoo = 1.7; a = 
Figure 20. Axial velocity contours 



(b) x = 1.032 in.; Uj/Uoo = 1.7; a = 4°. 
Figure 20. Concluded. 
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Figure 21. Vertical distribution of boattail axial velocity profiles for nonaxisymmetric boattail model, x = 0 in. 
y = 0 in. 
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Figure 22. Vertical distribution of boattail axial velocity profiles of nonaxisymmetric boattail model with solid 
plume included, x = 0 in.; y — 0 in. 
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Figure 23. Effect of angle of attack on boattail axial velocity vertical profile of nonaxisymmetric boattail. 
x — 0 in.; y = 0 in. 
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Figure 26. Reynolds normal stress (u'u'} profiles for nonaxisymmetric boattail model, y — 0 in.; a — 0 
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Figure 47. Streamwise vorticity contours. 
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